Environmentally-induced effects on a bipartite two-level system: geometric phase and 

entanglement properties 
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We calculate the geometric phase of a bipartite two-level system coupled to an external environ- 
ment. We analyze the reduced density matrix for an arbitrary initial state of the composite system 
and compute the correction to the unitary geometric phase through a kinematic approach. In all 
cases considered, we observe a similar structure as a function of the degree of the entanglement of 
the initial state. Further, we compute the entanglement entropy and concurrence of the bipartite 
state and analyze if there is any relation among these quantities and the geometric phase acquired 
during the nonunitary system's evolution. Finally, we discuss the results obtained. 
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I. INTRODUCTION 

It is plausible to imagine a quantum computer as a de- 
vice made of qubits which interact among themselves in 
some way. Therefore, a quantum computer can be under- 
stood as an interacting quantum many body system. In 
the last years, several results have been obtained which 
suggest that entanglement is at the root of the power of 
quantum computers [l[. 

A qubit is a quantum two-level system, i.e. a physical 
system described in terms of a Hubert space C 2 . This can 
be considered as a spin-1/2 particle, or as an atom with 
two energy levels or even a photon with two possible or- 
thogonal polarizations. Solid state examples which have 
achieved some experimental success include spins in semi- 
conductors and quantum dots, various designs based on 
superconductors, vacancy centers in diamond and singles 
molecule magnets (SMMs), see and references therein. 

It is has been shown that one can describe a quan- 
tum information processing system in terms of interact- 
ing qubits such as SMMs. Quite generally, in doing so, we 
need to include both environmental nonlocalized modes 
(like phonons and photons) and discrete localized modes 
(like defects, nuclear spins, loose spins). The spin-boson 
model corresponds to a single two-level system interact- 
ing with a large reservoir of bosonic field modes, i.e. a 
spin-1/2 particle coupled to an environment of harmonic 
oscillators. The seminal review paper by Legget et al. 
[H discusses its dynamics in great detail. The spin-spin 
model is its analogous but when the spin-1/2 particle is 
coupled to an environment of spin-1/2 particles. The ef- 
fect of this type of environment on the central spin has 
also been studied thoroughly by Zurek in Ref.Q. 
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When coherent superposition states on macroscopic 
scales are built, they are tipically created at very low 
temperatures in order to "freeze out" the thermal envi- 
ronment and thus minimize the decoherence effects. This 
thermal environment can be modelled into a bath of dclo- 
calized bosonic modes, implying an interaction among a 
large spatial region. But lowering the temperature does 
not affect the influence of localized modes such as nu- 
clear spins or impurities that are intrinsically present in 
the material Q. 

From another point of view, a system can retain the 
information of its motion when it undergoes a cyclic evo- 
lution, in the form of a geometric phase (GP), which 
was first put forward by Pancharatnam in optics [(| and 
later studied explicitly by Berry in a general quantal sys- 
tem 0- Since then, great progress has been achieved 
in this field. The geometric phase has been extended to 
the case of non-adiabatic evolutions Q. As an impor- 
tant evolvement, the application of the geometric phase 
has been proposed in many fields, such as the geomet- 
ric quantum computation. Due to its global properties, 
the geometric phase is propitious to construct fault tol- 
erant quantum gates. In this line of work, many phys- 
ical systems have been investigated to realize geometric 
quantum computation, such as NMR (Nuclear Magnetic 
Resonance), Joscphson junction, Ion trap and semicon- 
ductor quantum dots. The quantum computation scheme 
for the geometric phase has been proposed based on the 
Abelian or non-Abelian geometric phase, in which geo- 
metric phase has been shown to be robust against faults 
in the presence of some kind of external noise due to the 
geometric nature of Berry phase. It was therefore seen 
that interactions play an important role for the realiza- 
tion of some specific operations. Consequently, the study 
of the geometric phase was soon extended to open quan- 
tum systems. Following this idea, many authors have 
analyzed the correction to the geometric phase under the 
influence of an external environment using different ap- 
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preaches [9l4l5|. 

In this context, we shall briefly review the way the 
geometric phase can be computed for a system under the 
influence of external conditions such as an external bath. 
In Ref . Q , a quantum kinematic approach was proposed 
and the geometric phase (GP) for a mixed state under 
nonunitary evolution has been defined as 

<f> G = arg[ Vefc(0)£fc(r)(* fc (0)|* fc (r)) x 

^ k 

e - /J" !**,> J s (i) 

where are the eigenvalues and \^k) the eigenstates 
of the reduced density matrix p r (obtained after tracing 
over the reservoir degrees of freedom). In the last defini- 
tion, r denotes a time after the total system completes a 
cyclic evolution when it is isolated from the environment. 
Taking into account the effect of the environment, the 
system no longer undergoes a cyclic evolution. However, 
we shall consider a quasi cyclic path V : t e [0, r], with 
r = 2n/Q, (ft is the system's characteristic frequency). 
When the system is open, the original GP, i.e. the one 
that would have been obtained if the system had been 
closed <Pq, is modified. This means, in a general case, 
the phase is 4>g = 4>q + S<t>G, where 6(f>G depends on 
the kind of environment coupled to the main system 
Il2l . [HI UtT ] . If the eigenvalues of the density matrix are 
degenerate, the expression for the geometric phase takes 
a slightly different form, as described in ||. 

In this paper, we shall study the geometric phase ac- 
quired by a bipartite system in the presence of an ex- 
ternal environment. We shall consider both the presence 
of a bosonic and spin environment. We shall choose an 
arbitrary initial state and see how the geometric phase 
acquired during the evolution of the composite system 
is corrected by the presence of an environment. We will 
show the dependence of the environmentally induced cor- 
rection to the geometric phase upon the degree of entan- 
glement in the bipartite system. We will also show that, 
even in the nonunitary evolution, there is no correction 
to the phase when the composite system is in a maxi- 
mally entangled state (MES). In this case, we show that 
the system adopt a total phase <pQ = tt which, as it has 
been shown in previous articles, is of topological nature 
[l8j . Furthermore, we shall also study the entanglement 
measures of the bipartite system and analyze the relation 
among them and the geometric phase. 

This paper is organized as follows. In Section II, we 
study the evolution, geometric phase, and entanglement 
properties of the bipartite system in the presence of a 
bosonic environment at zero temperature. In Section 
III, we follow the same procedure for a spin environment 
leaving its trail in the dynamics of the bipartite system. 
Finally, we conclude our results in Section IV. Two ap- 
pendices complete the presentation. 



II. BIPARTITE SPIN-BOSON MODEL 

Oscillator environments correspond to a quasi contin- 
uum of delocalizcd bosonic field modes, with coherence 
and energy from the central system becoming effectively 
and irreversibly lost into this extended bosonic environ- 
ment. 

We shall consider a bipartite system, that is to say, two 
interacting two-level systems, both coupled to an exter- 
nal reservoir. 




Spin 1 




Spin 2 



FIG. 1: Diagram of the model we shall consider: two two-level 
systems coupled to an external reservoir whether it is bosonic 
or comprised of spins. 



In terms of the Hamiltonians, the model can be math- 
ematically described by the Hamiltonian of the free bi- 
partite system iJg, the Hamiltonian of interaction be- 
tween the bipartite and the external bath Hj and the 
free Hamiltonian of the external bath Hb ■ 

H s = -^crl + ^-at+jal®^ (2) 

N N 

Hi = <rl®^2 Kq n + cr 2 z ® ^2 g n q n (3) 

n— 1 n— 1 

N 

H B = ^2 hu n aia n , (4) 

n=l 

where the constants X n and g n couple the system to 
each oscillator in the environment, and 7 is the coupling 
strength between both spin-1/2 particles (as shown in 
FigOJ. Here we have assumed that each coupling con- 
stant of the two level systems with the environment is 
different being A„ for the spin I and g n for spin 2. 

In order to compute the geometric phase for the bi- 
partite system it is important to know its dynamics at 
all times. Therefore, in Appendix we have derivated 
the reduced density matrix of the bipartite system. It is 
important to know that this derivation has been done in 
the weak coupling limit for a general environment defined 
by a spectral function J (to). Following the mentioned 
derivation, if one consideres the most general case for an 
initial state of the bipartite system, namely 

|*(0)) = o|00) + P\0l) + C]10) + S\U), (5) 
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the reduced density matrix for this model can be written as, 
I 



Pr(t) 



( M 2 ap*e- i ^+^ t r 2 A 12 aCe-^+^^Au aS*e- l ^+ n ^ t T 1 T 2 T 2 12 \ 



^•e^+notrjiijj, c/?*e ?;(ni " n2) *rir 2 f 2 2 |c| 2 c<5*e-' i(02 " 27) *r 2 yi 



\6\< 



J 



where we have not written explicitly the dependence 
upon the time of Ti(t) just to simplify notation and we 
have defined 

f? 2 = f 2 2 (i) = c + 8 /o t *i Ji, »(*0 (6) 

and 

A12 =A 12 {t) =e i4 /o'««i G «(*0 (7) 

with Ti, r 2 , -Fi(t) and Gj(t) as defined in Appendix 
lAl Notice that the coupling constants of the model A ra 
and g n are absorbed in the dimensionless dissipative con- 
stants 701 ~ A 2 , 7o 2 ~ g 2 , and 7o 12 ~ Xg respectively, 
defined in the spectral density of the bath (sec Appendix 
[£]for details). 

In order to have a complete description of the effect 
of the bath on the bipartite system, we shall consider 
a quantum ohmic environment at zero temperature and 
different couplings between both spin-particles and the 
bath. Consequently, there will be different decoherence 
factors (see Appendix [SJ and we shall use subindcxes 
1,2, or 12 to refer to them. By the definition of the 
particular initial state, the reduced density matrix and 
consequently, the dynamics of the open bipartite system 
is known. 

For the particular purpose of this work, we shall de- 
fine an initial density matrix of the form (by setting the 
corresponding values of Eq. ([5]) ) 

ft(O) = ^/ + 7#)<0|, (8) 

where r e(0, 1] determines the mixing of the state and I 
is the unit matrix in the Hilbert space 2x2. The state 
\(j>) may be any of the following states, 

\&) = 0^100) + v$|ll) (9) 
= v^^|01) + Vp|10) (10) 

where p determines the degree of entanglement being 
|0), |1) eigenstates of the Pauli operator a z . It is easy to 



note that when p = 1/2, Eqs.© and (JTUJ) are Bell states 
and Eq.® defines the so-called Werner states which play 
an important role in quantum information processing. 
We use Eq.(j8]) because it includes all possible cases, such 
as pure or mixed states and maximal or non-maximal en- 
tangled states. The first term in Eq.® can be regarded 
as the noise and the mixing coefficient r describes the in- 
tensity of noise. Recently the one-to-one correspondence 
between r of the Werner state and the temperature T 
of the one-dimensional Heisenberg two-spin chain with a 
magnetic field B along the z-axis, has been established 
0. 

Finally, if we assume an ohmic environment at zero 
temperature, the decoherence factors take the following 
forms, 

rx(t) = e -270 1 log(l+A 2 t 2 ) j (u) 

r 2 (t) = e^WogU+A 2 ^ (12) 

r 12 (t) = e - 2 Toi 2 log(l+A 2 t 2 ) j ( 13 ) 

similarly to the spin-boson model for zero temperature 
[Icl ] , where A is the environmental frequency cutoff. 

At this stage, we know the dynamics of the bipartite 
system for all times. However, we are interested in 
the geometric phase adquired by the composite system 
in one quasicyclic evolution r ~ 2n/D,, being Q the 
characteristic frequency of the bipartite system. As 
we have mentioned in the Introduction, the kinematic 
approach to the geometric phase can be done by the 
use of the reduced density matrix. Therefore, we can 
compute the eigenvalues and eigenvectors of the matrix 
and compute the geometric phase of this state. 



A. Werner state with |#) 

We shall start by assuming an initial density matrix of 
the form of Eqs.([8]) and (|9]). In such a case, it is possible 
to write the reduced density matrix for all time t > as 
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In this case, the eigenvalues of Eq. (jT4")) are 
1 



e+(t) = ^(l + r + 2r v /l + 4p(l-p)(rf 2 r?rl-l)) J 



1 

4 

1 - r 

"I - ' 

1 - r 
~4~' 



e_(t) = |(l+r-2rVl + 4p(l-p)(rf 2 rfrl-l)), 



£l 
£"2 



with the corresponding eigenvectors, 
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+ P (i-p)r2rfr2rf2 



|00) 



11) 



|00) 



11) 



101) 
|10). 



(15) 



We see that, with this particular choice of the initial 
state, we are left to work in the space spanned by 1 00) 
and 1 11), since ^1,2) are time independent and hence 
will not contribute to the geometric phase. Using Eq.([T]), 
formally the geometric phase can be computed as 



4> G = arg| y/ e+ (0)e+ (r) <*+ (0) |*+ (r)) 

g i Q/ o T cos 2 (0 + (*i)) dti 



v/ e _(0) e _(r)(*_(0)|*_(r)) 
e i n/ T co B 2 (0_(ti)) dti 



(16) 



by defining 
cos(0±(i)) 



Vp(i-p)rir 2 r 



U/( e± -(i-p)) 2 +p(i-p)rf 2 r?ri)' 

(17) 

in analogy to what has been done for a single two-level 
system in SU(2) 0, [l(J. In this case, one can use 



r = 2ir/n, with O = (fii + fi 2 ). 

In order to compute the geometric phase and obtain 
concrete results, we shall choose r = 1. In this case it is 
easy to see that Eq. ([TB)) simplifies since £-(0) = and 
the contribution of this eigenenergy and its associated 
eigenvector to the geometric phase is null. Hence, the 
the geometric phase becomes 



(fii + fi 2 ) x 

dt 



(i-p)pr*,r?ri 



o (l~p)pTi 2 rm + [e 

(fii+fi 2 ) f cos 2 (0+(i)). 
Jo 



(18) 




FIG. 2: Deviation from the unitary geometric phase for the 
bipartite system coupled to an ohmic bosonic environment as 
a function of the initially entanglement of the state p and the 
coupling constant 70. The frequency cutoff is A/fi = 100. 

It is important to note that if the bipartite system is 
isolated, then Ti(t) = 1, and the geometric phase be- 
comes 

4> G = {tt 1 + Sl 2 ) ( dt (l-p) = 2vr(l-p), 
Jo 

which, of course, is the unitary geometric phase </>q. We 
can interpret this result by analogy to a single spin-boson 
model. In Ref.[l(3], we have started by a pure initial state 
on the Bloch Sphere, characterized by #0, i.e. |3>(0)) = 
cos(0o/2)|O)+sin(#o/2)|l). Then, in our present case, it is 
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cos(#o/2). Then the unitary 
>g = 27rcos 2 (6» /2) = tt(1 + 



enough to define ^/l — p ■■ 
geometric phase becomes 
cos#o), which agrees with the geometric phase acquired 
by a spin in SU(2). In other words, as the subspace 
spanned by {|01), |10)} is a decoherence free space, it is 
not affected by the presence of an environment. Since 
[Hs, Hi] = 0, in this case, it does not evolve in time (see 
Eq. (|Aip ) . Then, as long as the whole system is prepared 
in a separable state, the geometric phase of the composite 
is exactly that of the evolution happening in the Hilbcrt 
Space spanned by {|00), |H)}, similarly to that of 2 x 2. 

It is also important to note that in the case of having 
a maximally entangled state is 4>g — Physically, this 
phase depends on the parity of the number of times the 
state crosses the space orthogonal to the initial state in 
the representation of a S03 sphere. It is already known 
that for a qubit the total phase gained is 7r (or n it) and 
it is due to a combination of the dynamical phase and 
the geometric phase. For a two- level bipartite system 
with an arbitrary degree of entanglement a third possi- 
ble type of global phase can be identified (in addition to 
the dynamical and geometric phase already known): a 
topological phase, which is a consequence of the geom- 
etry of the entangled two-level system. This phase has 
been studied for maximally entangled states (MES) and 
it is at the origin of singularities appearing in the phase 
of MES during a cyclic evolution. In |l8|, it is studied 



the phase dynamics of entangled qubits under unitary 
cyclic evolutions. Therein, it is shown that, after a cyclic 
evolution, the combination of the different phases always 
leads to a global phase of an entire multiple of n. This 
result, already known and verified experimentally for a 
single qubit is recovered here for an entangled qubit with 
maximal degree of entanglement in the presence of an 
environment. 

For a bipartite state, we can not longer use the Bloch 
sphere to seek a geometric representation of that state. In 
[1 81 ] it has been shown that a geometric representantion 
of a bipartite state can be obtained by using a Bloch ball 
and a S03 sphere. Therefore, the total phase gained by a 
state is a combination of not only the dynamical and geo- 
metrical phase, but also the topological phase. Similarly 
to one qubit states, MES also gain a total phase of ir (or 
rm) under a cyclic evolution. However, this phase is of 
topological origin. In this context, we can explain our re- 
sults. We arc looking to corrections to the unitary phase. 
However, in the case of a MES there is no correction to 
the unitary phase. How can this be? We strongly believe 
that it is due to the nature of the unitary phase. As we 
have explained above, it is of topological origin. Then, it 
can not be disturbed or modified by the dynamic of the 
environment. For all other degree of entanglement, we 
do obtain a correction displaying the type of behaviours 
stated in [IH, with the interesting additive factor that 
we are considering the presence of an enviroment. 

In Fig|2j we show the behavior of the geometric phase 
for the open composite system, as a function of p and the 
coupling constant 70, when the bipartite system is cou- 



pled to a bosonic ohmic environment by the same dimcn- 
sionless coupling, i.e. 70 j = 7o 2 = 7012 = 7o- Therein, 
we can note that for very small values of the coupling con- 
stant, the geometric phase is that of the unitary system, 
meaning by the latter that there is no interaction of the 
bipartite with an environment (closed system). When 
the presence of the environment is relevant enough to in- 
fluence the dynamics of the bipartite system, the correc- 
tion to the unitary geometric phase 8<j)G also depends of 
the degree of entanglement p, being zero when the initial 
state is maximally entangled (p = 1/2). This result has 
so me p eculiar characteristics similar to the ones observed 
in [20| . There, authors studied the isolated system and 
obtained the geometric phase by tracing over one spin 
as suggested in plj . Even though we are considering an 
external environment, we can note the symmetry in p of 
our result derived from the type of coupling, and also the 
jump at the crossover point p — 1/2 (8q = 7r/2), when the 
initial state falls from the upper to the lower semi-sphere 
in the Bloch representation. 

In order to analitically analyze the correction to the 
unitary geometric phase we can perform a series expan- 
sion in powers of the coupling with the environment 70 
(assuming a solely coupling between the bipartite system 
and the environment, i.e. 70^ = 70)) 

2n(l -p) + 32 7o p^(l -3p + 2p 2 ) 



0G 



arctan(27r— ) 



A 

where Q, = f2i + O2 



los 



A 



2 A 2 



(19) 



Assuming that A/fi ^> 1, the correction 6(/>g to the 
geometric phase can be well approximated by 



( log (\-) 



(20) 



which is essentially the same correction found in fioj for 
a single two-level particle. 

In Figj3]we show the geometric phase from Eq. (fT8)) and 
the geometric phase obtained by the use of Eq. (fT9"|) . The 
above expressions show that the phase of the composite 
system depends on several parameters: the coupling 
constant 70, the frequency cutoff A (i.e. the "size" of 
the environment) and the degree of entanglement of 
the initial state p. In FigEl we can observe that the 
correction grows with p up to p = 1/2. For values of 
"fop <C 1 the perturbative expression for the geometric 
phase becomes a good approximation. 

We can think of a realistic model where our central 
system is coupled to an external environment composed 
of phonons or photons. In such case, we should con- 
sider an non-ohmic spectral density for such an environ- 
ment, particularly one that goes as J(ui) ~ w 3 , usually 
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A series expansion in powers of the coupling constant 
70, gives 



b G « 2^(l-p) + 8 7o p^(l -3p + 2p 2 ) 

TAI n + n» + 4^A0 " 3MCtaa(2T n ) 



0.01 

Yo 



0.02 



(22) 



where the correction to the phase can also approximated 
for A/O > 1 as 



FIG. 3: (colors online) Deviation from the unitary geometric 
phase for the bipartite system coupled to an ohmic bosonic 
environment as a function of the initially entanglement of the 
state p and the coupling constant 70. Exact results are plotted 
with lines and perturbative ones with lines and dots. Equal 
colors indicate equal values of p. A/fi = 100. 
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Scj) G « 327170^(1 - 3p + 2p 2 ), (23) 

showing it is smaller than the correction induced in the 
ohmic case Eq.® by a log(A/0) factor [ijj. 
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FIG. 4: (colors online) Deviation from the unitary geomet- 
ric phase for the bipartite system coupled to an supraohmic 
bosonic environment as a function of the initially entangle- 
ment of the state p and the coupling constant 70. A/f2 = 100. 

called supraohmic environment Q • The decoherence fac- 
tor takes a subtly different form, being decoherence even 
less effective at zero temperature, 

psupra^ = e - 4 7o (1+ A A 4 4^ _ ( 2 1) 

For times At > 1, this factor becomes r supra (t) ~ e~^°. 
Since 70 < 1, decoherence is not as effective as in previ- 
ous the cxa mple[l3[23. InFigH we show the geometric 
phase for the open bipartite system when it is coupled to 
a bosonic supraohmic environment by the same coupling 
constant 70. In this case, we see that the effect of the en- 
vironment over the geometric phase is less relevant since 
the correction to the unitary geometric phase is smaller 
than the previous case for the same set of parameters. 



FIG. 5: (colors online) Deviation from the unitary geometric 
phase for the bipartite system coupled to a nonohmic bosonic 
environment as a function of the initially entanglement of the 
state p and the coupling constant 70. Exact results are plotted 
with lines and perturbative ones with lines and dots. Equal 
colors indicate equal values of p. A/Q = 100. 

In FigJS] we present the deviation from the unitary 
geometric phase for the bipartite system coupled to a 
supraohmic bosonic environment as a function of the de- 
gree of entanglement of the initial state p and the cou- 
pling constant 70. The supraohmic environment is less 
effective in inducing decoherence on the system and con- 
sequently, affects less the system's dynamics, but it also 
produces a phase correction on the bipartite state. 

The correction to the unitary geometric phase natu- 
rally depends on the size of the environment, set by A. 
If we consider bigger or smaller ones, the result will be 
qualitatively the same [To|, [H|. Typically, decoherence 
from the ohmic bath is larger, limiting the possibility of 
measuring the correction to the phase using interferom- 
ctry. 
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1. Linear Entropy and Concurrence 

Entanglement is a quantum-mechanical feature which 
does not exist in the classical world. It carries non-local 
correlations between the different parts in such a way 
that can not be described classically. Bipartite entan- 
glement of pure states is conceptually well understood. 
A useful measure of many-body entanglement when the 
total system is in a pure state is the Von Neumann en- 
tanglement entropy Sl(t), which provides a measure of 
the bipartite entanglement present in pure states. To be 
precise, the entanglement entropy measures the optimal 
rate at which it is possible to distill Bell pairs by local 
operations in the limit of having an infinite number of 
copies of the bipartite system. 
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FIG. 6: (colors online) Concurrence and Linear Entropy as a 
function of time for different values of the parameter p when 
the bipartite system is coupled to a bosonic environment with 
coupling constant 70 = 0.002. Parameters used: A = 100, 
p — 0.01 solid line, p = 0.20 dashed lined and p — 0.5 dot- 
dashed line. The concurrence for a supraohmic environment 
(with dots) is also indicated for p — 0.01 and p — 0.5. 

The von Neumann entanglement entropy (or linear en- 
tropy) is obtained by focusing on bipartite systems where 
space can be divided into two regions: the one corre- 
sponding to the bipartite system and the one of the en- 
vironment. After tracing over the degrees of freedom of 
the environment, we obtain the reduced density matrix of 
the system p r . The von Neumann entanglement entropy 
is defined as, 

SI = -Tr[p T logM}. (24) 

In terms of the eigenvalues of the reduced density ma- 
trix of the bipartite system, the entanglement entropy SI 
reads, 

Sl(t) = - e+ (i)log 2 ( e+ (t)) -e_(t)log a (e_(i)). (25) 

Quantum dccohcrcncc implies a rapid reduction of the 
off-diagonal terms of the bipartite reduced density ma- 
trix which results in the case of maximal entanglement 



in e±(t) 1/2 and Sl(t) -> 1. 

The quantity for measuring the entanglement between 
the different parts of the composite system is the con- 
currence [23|. The concurrence for the evolution of this 
state can be computed as C(p Cr ) = max(0, \/Al — V%2 ~ 
V^3~ V'Ai), where Ai, A2, A3 and A4 are the eigenvalues 

In the case of the initial condition considered in the 
last Section, the concurrence is 

C A = 2r 4 (t)^/p(l- P ). (26) 

In the closed system case, the concurrence is Ca = 1 
when considering a MES (p = 1/2), and Ca = for 
product sates (p = 0). In the open system case, the con- 
currence depends on the decoherence factor T(t) as can 
be seen from Eq. (|26j) . being smaller in the supraohmic 
case. 

In Figs |6] and [7] we show the linear entropy (Eq. (|25|l ) 
and concurrence (Eq. (|26[) ) when the bipartite system is 
coupled to the bosonic environment by a solely coupling 
constant 70. Therein, we can see that when the value 
of 70 is small, then the ohmic and supraohmic environ- 
ments induce the same behaviour on the entanglement 
of the system state. However, as 70 increases, the ohmic 
environment destroys whatever degree of entanglement 
might be while the supraohmic is not effective in such a 
task. 

1 T- — , , , p=0 5 Q 

0.9 \ \ 

0.8 \\ \ / 




p=0.20 . 



p=0.01 . 



0.02 0.04 0.06 0.08 0.1 

fit 

FIG. 7: (colors online) Concurrence and Linear Entropy as a 
function of time for different values of the parameter p when 
the bipartite system is coupled to a bosonic environment with 
coupling constant 70 = 0.1. Parameters used: A = 100, p = 
0.01 solid line, p = 0.20 dashed lined and p = 0.5 dot-dashed 
line. The concurrence for a supraohmic environment is also 
indicated (with dots) for p — 0.01 and p = 0.5. 

It is worth mentioning that both quantities, i.e the ge- 
ometric phase and concurrence, are modified by the pres- 
ence of the environment to a greeter or lesser extend. We 
can see that for a strong environment (70 = 0.1), and en- 
tangled states between < p < 1/2 and 1/2 <p < 1, the 
correction to the geometric phase is important. However, 
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for all values of p the concurrence decays at very short 
times in the case of the ohmic environment. In the case of 
the supraohmic environment, the situation is much simi- 
lar to a quasi-isolated system, and the concurrence is not 
much affected by the noise and dissipation introduced by 
the external bath. 

Finally, it has been stated that there was some kind of 
correlation among the geometric phase of the individual 
spins and the concurrence for an entangled state of two 
spin- 1/2 particles [24[. However, we have shown that 
there is no such relations in the general case (when the 
coupling to an external environment is considered). Only 
in the isolated (or unitary) case, it is possible to explicitly 
show that 



C A 



l-p 



(27) 



which is constant for a fixed value of p. In the case of a 
MES with p = 1/2, the above relation reads $q = nC\ 

M. 



B. Werner state with \fi) 



Another possibility is to begin with an initial bipartite 
state of the form \y) in Eq.©. In this case, the reduced 
density matrix has a slightly simpler expression such as, 







(i- 





+ r(l -p) 



ry/p(l - p)e t ^-^ t T 1 T 2 f\ 2 






-p)e- < (" i -" 3 ) t r 1 r 2 f?2 o 

^+rp 

o fi? 



As before, we need to compute the eigenvalues and 
eigenvectors of the reduced density matrix in order to 
know the geometric phase of the system. For simplicity, 
we shall choose again r = 1. In that case, the eigenvalues 
will be very similar to those calculated before, with e\ = 
e 2 = and 

£+(*) = 5(1 + ^i + ^a-pXrVfri-i)), 
£-(*) = ^(1- Vi + ^i-pKf^rfri-i)), 

with the corresponding eigenvectors that can be easily 
computed as before. 

Again, we can see that £_(i = 0) = 0. Then, the 
only contribution to the geometric phase comes from the 
eigenvalue e + and its associated eigenvector |4 r +). Simi- 
larly to the above procedure, the geometric phase is, 

4> G = (Hi - 2 ) X 

[ T dt P (i-p)r;rlfj 2 
Jo p(i-p)r?rift 2 + [ e+ -(i-p)]a" 

From this equation, it is possible to show that there is 
no correction to the unitary phase in the case of starting 
whit a MES. Therefore, when p = 1/2 we get 4>g = 7r. 

As expected, we reobtain the unitary result for the 
geometric phase if the system is isolated, i.e. joi = 0, and 
consequently 1^ = 1. Furthermore, if we assume that the 
coupling between each spin and the environment is the 



same, i.e. 701 = 702 = I012 = 70i we see that surprisingly 
there is no correction to the unitary geometric phase, 
since 

cb G = 4> U G =2tt(1-p) 

being at this time r = 27r/Sl with fl = (Qi — SI2). 

We can state that in this case the geometric phase is ro- 
bust against the action of the external environment. This 
is easily understood since, as we have explained before, 
the subspace spanned by 1 01} , 1 10} is decoherence-free. 
Solely in the case in which all the coupling strengths are 
different, one to each other, the total geometric phase 
accounts for environmentally induced corrections. 

Incidentally, if we compute the concurrence for this 
initial bipartite system state, we see that it is no longer 
a function of the decohcrence factor (when 70^ = 70), 



C B = \jl-p + 2^p{l-pf 
- ^l-p-2^pJT~pf. 

(28) 

In this case the concurrence Cb does not depends on 
the decoherence factor and, as in the usual case, we have 
C B = 1 for p = 1/2 and C B = for p = 0. 
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III. BIPARTITE SPIN-SPIN MODEL 

Spin environments are typically the appropriate model 
in the low temperature regime. In particular, experi- 
ments devoted to the studies of macroscopic quantum 
coherence and dccohcrcncc require temperatures close to 
absolute zero in order to operate. Experimental evidence 
shows that dccohcrcncc is mainly dominated by interac- 
tions with localized modes in this setting, such as para- 
magnetic spins, electronic impurities, defects and nuclear 
spins. Each of the localized modes is described by a 
finite-dimensional Hilbert space with a finite energy cut- 
off. There are numerous examples that show how spin 
environments influence the dynamics of a central system 
by becoming a source of decoherence. 

In this framework, we examine the decoherence in- 
duced by disordered interacting spin baths at finite tem- 
perature. Our choice of the bath is the most simple case 
for an Ising chain so as to facilitate an analytical study of 
the decoherence and the geometric phase for the model. 
We shall consider the bipartite two-level system coupled 
to an external spin environment, modelled by the follow- 
ing Hamiltonians: 



We have included the free Hamiltonian Hb of the envi- 
ronment, where hi denotes the tunneling matrix element 
for the ith-environmcntal spin. This free Hamiltonian 
lends intrinsic dynamics to the environment, in contrast 
with the more simplified spin-environment models. 

As we have done with the bosonic environment, it is 
imperious to know the dynamics of the system in order 
to study the geometric phase from the kinematic point of 
view (Eq. (JTJ) ) . To that end, we have derived the reduced 
density matrix of the bipartite system in Appendix [B] 
Herein, we shall consider two different initial states as 
before (Eqs. © and (fTU|)), 



Hi 



-(n 1 a 1 z +n 2 a 2 z )+-,<j 1 z <Z>a 2 z 



N 



N 



i=l i=l 

N 



(29) 



A. Werner state con |#) 



Herein, we shall write in Eq.©. For this initial 
state, the reduced density matrix has a much simpler 
expression than the general one, 



Pr A (*) 








tVp(1 -p)e~ i ^ +n ^ t Q{t) \ 

4 











\ tVp(1 -p)e i( - n2+Ql ^Q*(t) 

I 



+rp 



where the decoherence factor Q(t) is defined as, 



Q(t) = n{i 

i=l ^ 



2(e l + X i y 



x sin 2 ^^ + (£,, + A;) 2 ) 



(30) 



and has been derived in the Appendix [B] 

In order to compute the geometric phase, we need to 
know the eigenvalues of the reduced density matrix. By 
setting r = 1, these are 



and the corresponding eigenvectors 

M = loi) 



and 



ei,2 = 
I 
2 



i±Vi + 4(i-pMQW 2 -i) 



(31) 



K) = 



W) = 1 10) 
Vp(l-p)IQ(*)|e _i(ni+n2)t 



Vp(i-p)\Q(t)\ 2 + 


[c+-(l- 


-p)f 


[6+-(l- 


-P)} 




^ P (i - P )\Q(t)\ 2 + 


[c+-(l- 


-p)f 


^p(l-p)\Q(t)\ 


e -i(a 1 +a 2 


)t 


^/p(i - P )\Q(t)\ 2 + 


[e- - (1 - 


-p)Y 


[e- (1 - 


-P)] 





^p(l ~p)\Q(t)\ 2 + [e--(l- P)f 



|00) 
|H) 
|00) 

111) 
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By the use of Eq.([T|), the geometric phase becomes, 

p(i- P )\Q(t)\ 2 



6 G = (fi 1 + 2 ) 



dt 



p(l-p)\Q(t)\2 + ie + -(l 



P)} 2 ' 
(32) 

and it is represented in Fig. (jHJ for an environment of 
A" = 100 spins. Therein, we can see that the deviation 
from the unitary geometric phase is less than in the ohmic 
bosonic environment seen before, but greater than in the 
supraohmic case. Then, if we consider an experiment of 
a bipartite system coupled to a spin and photonic en- 
vironment, we can note that decoherence comes mainly 
from the spin environment as state in Q (assuming "size- 
comparable" environments) . 




dimensionlcss coupling constant X/h. In such a case, 
4> G « 27r(l -p) 

2 



16ATp(l - 3p + 2p 2 



n ■ 




47T — - sin 


SO] 


h 





(34) 



where f2 = fii + Q 2 
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FIG. 9: (colors online) Deviation from the unitary geomet- 
ric phase for the bipartite system coupled to an environment 
comprised of N=100 spins as a function of the entanglement 
of the initial state p and the coupling constant X/h. Exact re- 
sults are plotted with lines while perturbative ones with lines 
and dots. Equal colors indicate equal values of p. 



FIG. 8: (colors online) Geometric phase for a bipartite system 
coupled to a spin environment comprised of N = 100 spins. 

As we have shown for the bosonic environment, the 
correction to the geometric phase is zero for an initial 
MES of the composite system, getting a total phase of w. 

It is easy to note that we can reobtain the unitary 
geometric phase if the bipartite system is isolated. In 
that case, e, = = A^ in Q(t) and the decoherence factor 
Q(t) = 1. Similarly, we can check that e+ = 1 and 
e_ = 0. Then the geometric phase becomes 

G = (Ox + 2 ) f dt(l -p) = 4> U G = 2rr(l - P)> ( 33 ) 
Jo 

assuming r = 27r/(f2x + ^2)- 

We can perform a series expansion in powers of the cou- 
pling with the environment. For that, we shall assume 
that the couplings between each spin and the environ- 
ment are equal, i.e. £{ = X{. In order to achieve an ana- 
lytical result, we shall assume all the bath spins to have 
the same coupling constant and frequency Then, the de- 
coherence factor becomes Q(t) = Y\f =1 q%{t) = q{t) N and 
we can forget about the product function. Besides, we 
shall perform a perturbative expansion in powers of the 



In Fig|9] we can see that the series expansion of the 
geometric phase results a good approximation of the 
latter for mostly all values of p. It is more stressed in 
the case of small pN since the correction to the phase is 
quadratic in X/h. 



1. Linear Entropy and Concurrence 

As in the other examples analysed for the bosonic envi- 
ronment, we can study the behaviour of the concurrence 
during the dynamics of the bipartite system. In this case, 
if we compute the concurrence as explained above, we ob- 
tain 

C A = y/p(l - P )(Q(t) + l) 2 - y/p(l - P ){Q{t) - l) 2 . 

(35) 

In FigflQl and HH we plot the concurrence and linear 
entropy as a function of the time for different values of 
p and the coupling constant A = e, for an environment 
comprised of A^ = 10 and N = 100 spins. In both Fig- 
ures, we can see a peculiar behavior for some of these 
values. As expected, since we are dealing with a finite 
environment, information can be in principle recoverable 
since Q(t) is at worst quasiperiodic (both proportional 
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FIG. 10: (colors online) Concurrence as a function of time for 
different values of the coupling constant A and different values 
of the initial state p, when the bipartite system is coupled to 
an environment comprised of N = 10 spins. The dotted line 
represents the concurrence for p = 0.01, A = 0.1 and N = 100. 
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FIG. 11: (colors online) Concurrence and Linear Entropy as 
a function of time for A = 0.1 and different values of the 
initial state p, when the bipartite system is coupled to an 
environment comprissed of N = 10 spins in the case of p = 
0.45. For p = 0.1, we present both: an environemt comprised 
of N = 10 spins (red solid line) and N = 100 spins (grey 
dotted line). 



to the decoherence factor Q(t)). As it is well known for 
this type of environment, the effectiveness of the deco- 
herence mechanism is determined by the dimension of 
the environment. That is why we observe the periodic 
behaviour in the analyzed quantities. On the other size, 
when the coupling with the environment is very small, 
its presence is not of great importance in the dynamics 
of the bipartite, and the concurrence is then constant. 
However, as we increase the coupling constant A, we find 
that the concurrence oscillates due to the oscillating func- 
tion present in the decoherence factor. In Fig[TU]we have 
plotted the concurrence for different values of p (degree of 
entanglement of the initial state) with the same color for 
equal values of the coupling constant A. All values consid- 
ered with lines are for an environment of N = 10 spins, 
while the dots represent the concurrence for p = 0.01 
and A = 0.1 of an environment of N = 100 spins. We 
can observe that the effect is stronger when the environ- 
ment is bigger. Just for the sake of completeness, we 
plot in Fig llll the concurrence and linear entropy for two 
extreme values of the degree of entanglement p for a cou- 
pling constant A = 0.1. In the case of p = 0.1 we can 
again observe the difference in both quantities when the 
size of the environment is N = 10 (solid line) and when 
it is N = 100 (grey dotted line). 



B. Werner state with 



Herein, we shall write in Eq.®. For this initial 
state, the reduced density matrix has a much simpler 
expression, 



PrB(i) 







(1-r) 



"(1-P) 









rVp(l -p)e" i(f2l+f22)t P(t) 

i + rp 

^ 



u- 



where the decoherence factor P(t) is has been defined in the Appendix IBI 



N 



P(t) — Y\. \ 1 — ( }j2 \ f £ . _ \.)2 ) I n order to compute the geometric phase, we need 

to know the eigenvalues of the reduced density matrix. 



i=l 

x sin^i^ + fo-Ai) 2 )}. (36) 
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These are (for r = 1) 

Cl.2 = 



f ± 



1± Vl + 4p(l-p)(|P(t)| 2 -l) 



(37) 



and the corresponding eigenvectors that again can be eas- 
ily computed. 

By the use of Eq.([T]), the geometric phase becomes, 



n 



dt 



p(l-p)\P(t)\< 



o p{\ -p)\P(t)\* + [e+-(l- p)] 2 ' 



where SI = — Sl 2 ). It is easy to note that we can 
reobtain the unitary geometric phase if the bipartite sys- 
tem is isolated or if it is in a MES with p = 1/2. In 
isolated from the environment case, £, = = Xi in P{t) 
and the decoherence factor P(t) — 1. Similarly, we can 
check that e + = 1 and e_ = 0. Then the geometric phase 
becomes 



G = fi / di(l -p) = ^ = 2tt(1 - p), (38) 
Jo 

assuming r = 27r/0, once again. 

The concurrence for this case can be also calculated as 



c B = vWi - + 1) 2 - vW - - 1) 2 - 

As it is easy to note, this case is similar to the one 
presented before but changing the decoherence factor 
Q(t) —> P(t). However this fact is not trivial. If we 
consider that the coupling between each spin of the 
system and the external reservoir is equal, i.e. £j = Aj, 
the P(t) = 1 and the subspace generated by the states 
1 01) and 1 10) becomes a decoherence free subspace, reob- 
taining t he unita ry geometric phase and a concurrence 

C B = 2y/ P (l- P ). 



IV. CONCLUSIONS 

The geometric phase of entangled states is an issue 
worth of attention. It could be a potential application in 
holonomic quantum computation since the study of en- 
tangled spin systems effectively allows us to contemplate 
the design of a solid state quantum computer. However, 
decoherence is the main obstacle to overcome. All real- 
istic quantum systems are coupled to their surroundings 
to a greater or lesser extent. 

We have thoroughly studied the geometric phase for 
a bipartite system, i.e. two coupled qubits, also cou- 
pled to an external environment. We have considered 
both cases: a bosonic environment at zero temperature 
(whether ohmic or supraohmic) and a spin environment. 
In all cases, we have chosen a general initial state and 
computed the evolution of the composite system in terms 
of the reduced density matrix. We have further estimated 



the decoherence factors for all cases in order to have a 
full insight into the decoherence process induced by the 
environment on the system's dynamics. We have seen 
that there is a hierarchy among the environments. That 
means, when considering real systems, such as single 
molecule magnets (for example Fes molecule), decoher- 
ence induced by the nonlocalizcd modes of a supraohmic 
environment is less than that induced by the localized 
dicrctc modes of a spin environment. 

By the use of a kinematic approach, we have computed 
the geometric phase 4>o = 4>g + &4>g f° r different choices 
of the initial state of the bipartite system. In some cases, 
we have also performed a perturbative expansion of the 
geometric phase in order to complete the analysis. In 
all cases, we have found the same geometric dependence 
upon the initial angle 9q for the correction of the geomet- 
ric phase that has been found by some authors by the use 
of different approaches [l3j . 

As entanglement is considered to be one of the key 
resources in quantum information science, we have also 
studied this property of the composite system. In partic- 
ular, we have checked that the correction to the geometric 
phase is 54>g = f° r the case of maximally entangled ini- 
tial states of the bipartite systems. This is the case for all 
the cases considered, no matter the kind of environment 
would be present so far. 

We have also computed the concurrence and linear en- 
tropy of each initial state. We find a steady relation be- 
tween the geometric phase and the concurrence only for 
a Bell state in the isolated situation. In all other cases, 
not only the geometric phase but also the concurrence 
are modified but the presence of the environment. 

In all cases it is posible to compute the geometric phase 
in not only one cyclic but many, i.e. r = n27r/f2, with 
n integer. In that case, the correction to the unitary ge- 
ometric phase would be proportional to n, the winding 
number 12, EH- This proportionality of Sfic to n might 
imply that the geometric phase of the composite system 
still has some of its geometric character. However, this 
correction is also a function of the environments spec- 
trum and then the total geometric phase is not a simple 
geometric quantity. 



Appendix A: Derivation of the reduced density 
matrix for a bosonic environment 



Herein, we shall derive the reduced density matrix for 
a two spin 1 /2 particles coupled to a bosonic environment 
through coupling constants X n for the spin 1 and g n for 
spin 2 (being A„ not necessarily equal to g n ). 

In order to know the dynamics of the bipartite sys- 
tem, we shall find the reduced density matrix for the 
composite system, as has been done in [25| for one par- 
ticle spin-boson model in the case of weak coupling. We 
must note that in this case the interaction is V(ii) = 
a l En \iQn{h) + a 1 J2 n 9nQn(ti)- Similarly to the sin- 
gle spin-boson model, it is straightforward to compute 
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the time derivative of the reduced density matrix as (we 
have set H = 1 from here up to what follows) : 



dtiTrB 



Vfti), 



V(ti),/5r®p B (0) 



The tilde indicates that we are working in the Interaction 
Picture. 

After doing some algebra, the master equation for the 
reduced density matrix can be written as 

Pr(t) = -i[H s ,p r ] 



dh\ vi(ti)[<xl, [al,p r ] + v 2 {ti)[al, [o*,p r ]] 



+ ([(rl, Wlpr}} + [al [<rl,Pr]]) | 

+ ij dti Vi2(ti)([o-l,{(7l,p r }] 

+ WlWlpr}}), (Al) 



In Eq. (|Al[) we note the first difference with the de- 
phasing single spin-boson [l(J H^|. In the case of only 
one spin coupled to an external environment, there is 
no dissipation on the main system induced by the envi- 
ronment (it is just a pure dephasing model). However, 
in the bipartite system we see that dissipation appears 
dispitc of the similar coupling between the system and 
the environment. 

In this context, we can define the decoherence factors 
I\(f) as 



r i (t)=e- 4 /o«*i*'«(*0 ) (A2) 



with 



Fi(t) = / d£vi{t') 
Jo 

and the dissipation induced by the environment as 



G 12 (t) = f dt' m2 (t'). 
Jo 



where we have already taken the continuum limit and 
defined the noise and dissipation kernels as: 



Appendix B: Derivation of the reduced density 
matrix for an environment comprised of spins 



Vi{t) = / dui Ji{uj) cos(wi)coth( — ), 
Jo 2 



T)i(t) 



d(jjJi(uj) sin(ajt), 



with i = 1,2 or i = 12. Ji(u>) is the spectral density 
of the environment associated to each spin of the sys- 
tem, and Vi (t) the corresponding noise kernel while r\i (t) 
is the dissipation kernel associated to Ji{uj). One as- 
sumption we shall make is that J(lu) is a reasonably 
smooth function of oj, and that is of the form of 1 up 
to some frequency A that may be large compared to f2i 
and VL 2 . The spectral density function can be written 
as Ji{w) = 7 0l /4 w"A"- 1 e- w / A 0. Notice that at this 
stage, the coupling constants of the model A„ and g n 
have been absorbed in the continuous limit and that in- 
formation is now contained in the dimensionless dissipa- 
tive constants 701, 7021 an d 7012 respectively, defined in 
the spectral density of the bath. 



Herein, we shall derive the reduced density matrix for 
a bipartite system coupled to an environment comprised 
of N spins. The system's, environment's and interac- 
tion hamiltonians have been considered in the main text, 

secnni 

Let's take an initial separable state of the complete 
system, i.e. |*(0)) = |$(0)) <g> |x(0)), where |$(0)) is the 
initial state of the bipartite system while |x(0)) is for the 
spin environment. We can note that [Hs, Hi] = 0, so we 
can write the evolution of the complete state as 



(Bl) 



-i(H B )t 



x(o)>. 



After cumbersome calculations (that we do not con- 
sider necessary to be presented here), we have obtained 
the following reduced density matrix for a general initial 
state |$(0)) = a|00)+/3|01)+C|10)+£|ll) of the bipartite 
system 



/ \a\ 2 Q! ^ e -i{o a +2 7 )i M ^ 



^* e -i(n 1 -o 2 )tp^j ^5* e -i(^x-^)tR{t) 

ICI 2 C^e-^-^Sit) 
SCe^-^S^t) \S\ 2 J 



In the above expression, the trace over the degrees of freedom of the spin environment is contained in the 
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decoherence factors M(t), N(t), P(t), Q(i), R(t), 
and S(t). These factors have several expres- 
sions. M(t) = Tr B [e-^ HB + v ^ t pB(0)e i< - HB+v ^% 
with V a = (ei + Aj) and V r = (e,; — AA Sim- 
ilarly, N(t) = Tr B [e- l (^ +v '=)Vs(0)e 4(ffB - v ''- ) *], 
S(t) = TT B [e-^ HB+v ^ t PB {0)e^ HB - v ^ t ], 

R(t) = Ti B [e^ HB + v ^ t p B (0)e^ HB - v ^ t }, 

Q{t) TT B [e~^ HB+v ^ t p B (0)e^ HB - v ^ t ] and 

P(t) = Tr B [e~ l ( HB + v "*> t p B (0)e l ( HB - v ^ t }. In this 
paper, we shall only be interested in the last two coeffi- 
cients. Since we can choose a pure state for the initial 
state of the bath as |x(0)) = n^=i( a i|0«) + ft|li))> then 
Pb(0) = lx(0))(x(0)|- In order to have closed expressions 
for these coefficients, we have finally assumed \oti\ = \f3i\ 
and obtained the following decoherence factors 

x sin^i^f + ^ + Ai) 2 )} 
x sin^i^f + ^-Ai) 2 )}. 



We can find some common features between these ex- 
pressions and preexisting Literature. For example, we 
can note that the factor P(t) that affects the states |01) 
and 1 10) is similar to that found in [ljj for a two-level 
system coupled to an external environment. In that same 
case, if we do not consider the self interaction of the bath, 
i.e. hi = 0, then we get P(t) = Jl^li cos[2(ei + Aj)i] 
as found in 0, [H|. It is also possible to check that 
P(0) = 1. As it was done in we can evaluate 
the mean value of the decoherence factor. As in that 
case, (P(t))T^oo ~> 0. If we estimate the average dis- 
persion as a 2 = (P(t) 2 ) - (P(t)) 2 = E'iLiPh with 
Pi = (1 — jn£^^yX2\ )■ Under the assumption that 
all pi are approximately equal, the average fluctuations 
from zero are 

(B2) 

Therefore, large environments can effectively induce de- 
coherence on the central spin system. We can do the 
same with all the decoherence factors but in the follow- 
ing we shall only be interested in factors P(t) and Q(t). 
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